First-principles study of cubic perovskites SrMC>3 (M = Ti, V, Zr and Nb) 
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Using the full-potential linearized-augmented-plane-wave (FLAPW) method, we have analyzed 
systematically the trends in the structural and electronic properties of the 3d and 4d transition- 
metal oxides SrMOa (M — Ti, V, Zr and Nb). The optimized lattice parameters, bulk modules, 
densities of states, band structures and charge density distributions are obtained and compared with 
the available theoretical and experimental data. The energy gap between 02p - Md bands increases 
as the covalency of the system decreases going from 3d to 4d based perovskites. The electron 
configurations of Sr(Ti,Zr)C>3 and Sr(V,Nb)03 usually referred to as d° and d 1 oxides, respectively, 
differ considerably from these idealized "ionic" configurations, and the deviations increase with 
increasing of the d-p covalent overlap in the oxides. 
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I. INTRODUCTION 

Perovskites SrM03 attract much attention from 
physicists and material scientists because of the unusual 
combination of their magnetic, electronic and transport 
propertiesi2iiiiSi&. However, the electronic structure 
features in these complicated systems are not under- 
stood comprehensively yet, and attempts to explain 
their behavior in different conditions are often based 
on oversimplified models. Through quite a number of 
publications are known on strontium oxide perovskites 
with 3d metals on M siteaW^VMVWVe,! W9, 
the research focus has long been on strontium titanate, 
SrTiOs - one of the generic representatives of 3d 
transition-metal oxide perovskites. For instance, the 
LCAO-based approaches are used in non-self-consistcnt 
or self-consistent calculationsi^iSiii. Recently more 
accurate self-consistent full-potential LAPW, LMTO, 
ab-initio pseudopotential band structure calcula- 
tions are performedS&SiiSii^iiSiiS. Some attention is 
given to the related materials, viz., the 3d 1 metal- 
like perovskite SrV03, solid solutions SrTii-^V^Os, 
SrVOs/SrTiOa superlattices and heteroepitaxial struc- 
tures SrVOa/SrTiOa/SiSS^. Also, the cubic SrV0 3 
is treated as a model system for discussion of orbital 
fluctuations^ and correlation effects^. As an example of 
4d based perovskites, the cubic SrZrOs is studied by the 
discrete variation (DV) molecular orbital method 24 and, 
recently the band structure calculations are performed 11 
within the local-density approximation (LDA). 
Finally, the electronic properties of the ideal cubic 
SrNb03, which may be considered as a prototype of an 
ample family of perovskite-related layered niobates 2 ^, 
are investigated in workSS by means of LMTO-ASA 
approach. In this Communication, we report results of a 
systematic study of electronic properties of perovskite- 
type 3d and 4d transition-metal oxides SrM03 (M = Ti, 
V, Zr and Nb) by means of the FLAPW method within 
the generalized gradient approximation (GGA). We 
have evaluated a number of physical parameters such as 
optimized crystal structure, bulk modulus, density of 



states, band structure and electron density distribution 
and considered their variations with M-type metal. 

II. MODELS AND METHOD 

The considered SrM03 perovskites are assumed to 
have ideal cubic structure (s.g. Pm3m) where atomic 
positions in the elementary cell are M: la (0,0,0); O: 
3d (0,0,1/2); and Sr: lb (1/2,1/2,1/2). The electronic 
configurations are taken Ar3d 2 4s 2 for Ti, Ar3d 3 4s 2 
for V, Kr4d 2 5s 2 for Zr, Kr4d 3 5s 2 for Nb, Kr5s 2 for 
Sr and He2s 2 2p 4 for O. Here, the noble gas cores are 
distinguished from the sub-shells of valence electrons. 
The equilibrium structural parameters and bulk modules 
are calculated using the Vienna package WIEN2k, which 
is an implementation of the hybrid full potential (linear) 
augmented plane-wave plus local orbitals (L)APW+lo 
method within the density-functional theory 2 ''. The 
basis set inside each MT sphere is split into core and 
valence subsets. The core states are treated within the 
spherical part of the potential only, and are assumed to 
have a spherically symmetric charge density confined 
within MT spheres. The valence part is treated with 
the potential expanded into spherical harmonics up to 
1 = 4. The valence wave functions inside the spheres 
are expanded up to 1 = 12. The plane- wave expansion 
with RMT*KMAX equal to 10 and k sampling with 
10*10*10 k-points mesh in the Brillouin zone (BZ) is 
used. All calculations are carried out with optimized 
lattice constants. We use the Perdew-Burke-Ernzcrhof 
generalized gradient approximation (GGA)2 8 for the 
exchange correlation potentials. Relativistic effects 
are taken into account within the scalar-relativistic 
approximation. 

The self-consistent calculations were considered to have 
converged when the difference in the total energy of 
the crystal did not exceed 0.1 mRy as calculated at 
consecutive steps. We have adopted the MT radii of 1.7 
a.u. for O, 1.8 a.u. for M, and 2.5 a.u. for Sr. The 
lattice constants and bulk modules are calculated by 
fitting the total energy versus volume according to the 
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TABLE I: Calculated and experimental values for the lat- 
tice parameters, bulk modules and its pressure derivative of 
SrM0 3 (M= Ti, V, Zr and Nb). 



Parameters 


Lattice 
lattice (A) 


Bulk modulus 
B (CPA) 


Bo' 




Calcul. 


Exper. 


Calcul. 


Exper. 




SrTiOa 


3.9456* 
3.878 11 


3.90E 32 


167.9* 
191*1 


183^ 


4.56* 


S1-VO3 


3.8662* 


3.84l2i 
3.842C 33 


182.8* 




5.11* 


SrZrOs 


4.1794* 
4. 095 11 


4.109^ 


158.28* 
17lH 


150^ 


3.66* 


SrNbOa 


4.0730* 




170.0* 




3.71* 



* - our FLAPW-GGA data 



Murnaghan's equation of stated. The total density of 
states (DOS) was obtained using a modified tetrahedron 
method^. 

III. RESULTS AND DISCUSSION 

The lattice constants of the selected perovskites 
SrM0 3 , where M = Ti, Zr, V and Nb, were calculated 
under constrain of cubic symmetry. The results and 
some available experimental and theoretical data are 
summarized in Table 1. It is seen that for SrTiOa and 
SrZrOs the calculated LDA constants^ are smaller 
compared to experimental ones. The GGA improves the 
lattice constants over the LDA. In fact, the obtained 
GGA lattice constants are very close to experimental 
data; the deviations for SrTiOa, SrZrOs and SLVO3 
achieve only about 1%. The bulk modules decrease going 
from Si'Ti0 3 to SrZr0 3 and from SrV0 3 to SrNb0 3 ; and 
Bo values for metal-like SrV03 and SrNb03 are higher 
than for insulating SrTiOs and SrZrOs perovskites. 
These trends correlate with changes in the lattice 
constants. 

The dispersion E(k) in the high-symmetry directions 
in the Brillouin zone are shown in Fig. 1 where similarity 
in the energy bands is evident in the oxide perovskites 
under study. The general features in the spectra can 
be specified in more detail with the using of SrTi03 
as example. The lower bands, which are not shown in 
Fig.l, contain 12 electrons in Sr4p and 02s semi-core 
states located about -15.5 eV below the Fermi level. 



The upper valence band (VB) with the width of about 
4.6 eV is derived mainly from 02p orbitals with some 
admixture of Ti and Sr states, and it contains 18 
electrons. The bottom of the conduction band (CB) 
is composed mainly of Ti3d t2 g and e g states, and at 
higher energies the Sr4d - like bands are placed. The 
valence-band top of SrTiOs is composed of 02p orbitals, 
and it is located at the R point, while the VB maximum 
at the M point is about 0.09 eV lower than the R. The 
conduction-band bottom is at the T point. The CB 
minimum at X is by 0.11 eV higher than at the T point. 
Thus the indirect gap between R and T points is 1.81 eV, 
and the direct gap at T is 2.16 eV. These features of the 
band structure agree well with the previous calculations, 
see Table 2. At the same time the experimentally 
determined indirect band gap in SrTiOs achieves about 
3.30 eV, while the direct band gap energy is equal to 
3.75 eV^21. The divergence is related to the well-known 
underestimation of the band gap values within LDA - 
based calculation methods that treat electron exchange 
correlations in an approximate way. Notice, however, 
that the difference in the calculated values for the direct 
and indirect gaps equals exactly to the same difference 
in the experimentally observed band gaps. 
The nine occupied bands at the T point consist of three 
triply degenerate levels (Fis, T25 and Tis) with energies 
of -2.82, -1.17 and -0.39 eV below E F , respectively. The 
splitting of the levels (Fi 5 - r 2 5) and (r 2 5 - Ti 5 ) is due 
to the crystal field effects and electrostatic interaction 
between 02p orbitals. The triply (V 25 , +1.81 eV) and 
the doubly (Fi2, +4.12 eV) degenerate levels in the 
conduction band originate from t 2g (tt*) and e g (a*) 
bands of titanium, and they are separated with the 
energy gap of 2.31 eV. 

The calculated total (TDOS) and partial (PDOS) 
densities of states are shown in Fig. 2 within the energy 
interval [E F - 6 eV + B F + 8 eV]. The TDOS is almost 
entirely composed from the PDOSs for Ti3d and 02p 
states. The valence TDOS originates predominantly 
from 02p states, which is consistent with SrTiOs being 
a formal d° system with titanium in 4+ oxidation 
state. There is also a small contribution in the valence 
TDOS from strontium 5s states. However, there is 
a Ti3d PDOS in the occupied energy range due to 
the hybridization with the 02p states indicating the 
presence of a covalent bonding. Note also that within 
the cubic symmetry of SrTiOs the hybridization between 
t 2g and e g states is forbidden, and the orbitals within 
both bands are degenerate. Nevertheless, because of 
the Ti-0 hybridization the admixture of the both t 2 g 
and e g states to the oxygen 2p bands and also their 
partial overlapping occur. It is seen from Fig. 2 that the 
eg states (directed towards oxygen atoms so that they 
form a bonds with respective O 2p orbitals) have lower 
energy in comparison with Ti t2g states forming 7r-like 
bonds with O 2p states. 

The TDOS profile is formed basically of four peaks 
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TABLE II: Calculated and experimental values for the band- 
widths (in eV) of SrM0 3 (M= Ti, V, Zr and Nb). 



Parameters 


SrTiO 3 


SrZrOa 


SrV0 3 


SrNbOa 


Valence band 


4.59* 
4.9E 11 


3.86* 
4.32 11 


7.18* 
7.512 


8.12* 
6.132a 


Conduction band 
(occupied part) 






1.05* 
l.l^ 


1.21* 

1.182a 


Band gap: 
direct(r) 

indirect (R-r) 


2.16* 
2.3C 11 
3.7E 36 
1.81* 
1.9311 
3.2E 36 
3.421 


3.75* 
3.62±1 

3.30* 
3.3?li 
5.621 


1.40* 
1.02* 


2.79* 
2.38 

2.312a 



* - our FLAPW-GGA data 

* - Experimental data^2i 



centered at - 3.1, -1.3, + 3.6 and + 5.1 eV. A very 
flat 02p - like band (along T-X-M-r, see Fig. 1) is 
responsible for the highest of them placed at - 3.1 eV. 
The contribution of Ti3d states is quite small in the 
peak near -1.3 eV. Therefore the DOS peak at this 
energy may be attributed to quasi-flat bands originating 
primarily from non-bonding 02p states. The peaks 
in CB centered at + 3.6 and + 5.1 eV are formed by 
significant contributions from t 2s and e g bands of Ti3d 
states. 

The comparison of the band structures of the isoelec- 
tronic and isostructural SrTiC>3 and SrZr03 shows that 
the dispersion of the bands in SrTiC>3 is much stronger 
than in SrZr03. For example, the VB bandwidth in 
SrTiC>3 is about 0.73 eV larger. This means that Ti-0 
bonds are more covalent than respective bonds between 
zirconium and oxygen. Another obvious difference is 
that the band gap (which remains indirect in SrZrC^; 
the bottom of the conduction band is located at the T 
point in the cubic BZ while the top of the valence band 
is at the R point) is about 1.49 eV (45%) larger than 
that of SrTiC>3. In the case of SrZrC>3 the splitting of 
the bands (ris - T25) and (r25 - Tis) is about 0.75 and 
0.62 eV, respectively. 

Band structures for d 1 perovskites SrVC>3 and SrNbC>3 
are presented in Fig.l, and band structure parameters 
are listed in Table 2. The VBs of SrV0 3 and SrNb0 3 
consist of completely occupied oxygen 2p bands and 
partially occupied (V,Nb)3d bands. The electronic 
states near the Fermi level originate mainly from 3d 
(t2g) states. 



TABLE III: Occupation indices (in e) of M e 9 and t2 9 bands 
of SrMOs (M= Ti, V, Zr and Nb). 



Band/system 


SrTiO 3 


SrZrO a 


SrV0 3 


SrNbOa 


M e 9 


0.642* 


0.742 


0.307 


0.533 


M t 2s 


0.473 


1.313 


0.234 


0.791 



One can see from Fig. 2 that for SrV03 and SrNb03 
there is some contribution of the t 2g and e g bands in the 
energy intervals from -7.2 to -2.1 eV and from -8.3 to 3.6 
eV, respectively, because of the hybridization with 02p 
states. In the case of SrV0 3 the splitting of the bands 
(ri 5 - r 25 ) and (r 25 - ris) are about 1.94 and 0.81 eV; 
for SrNbOs these values are about 0.95 and 0.74 eV, 
i.e. larger than for SrTiOs and SrZr03, respectively. 
For SrV03 and SrNbOs the t2 S states are dominant in 
the vicinity of the Fermi level. The e g and t 2ff bands 
are well separated in the CB, and the maximum of the 
DOS in the eg band is centered above the upper band 
edge of the t 2ff band. The main weight of t 2g and e g 
states in SrV03 is observed in the energy intervals [1.1 
-i- +1.41 eV] and [+1.38 -f- +5.8 eV], respectively, while 
in SrNb03 the respective regions are located at [-1.2 ~ 
+2.3 eV] and [+1.61 ~ + 8.5 eV]. 

The obtained results show that the total density of 
states at the Fermi level N tot (E F ) in SrV0 3 (1.678 
states/eV/cell) is about 20% higher than in SrNbOs 
(1.256 states/eV/cell). The main contribution to 
Ntot(Ej?) in both systems originates from the t 2 g bands, 
and the N tot (Ej?) changes are also controlled by the 
contributions from t 2g states: going from SrV03 to 
SrNb03, these contributions decrease rapidly from 
1.109 to 0.542 states/eV/cell. In order to illustrate 
the bonding picture in systems under consideration 
we have plotted the charge density maps for SrTi03 
and SrV03, Fig. 3. It is seen that the covalent inter- 
action between Ti and O or V and O atoms is very strong. 

Finally, the perovskites SrTi03 and SrV03 (and 
their 4d analogues) are usually referred to as d° and 
d 1 oxides in the literature. It is seen, however, from 
Table 3 that calculated electron configurations (viz., 
Ti3d 1115 , Zr4d 0541 , V3d 2 055 and Nb4d 1 324 ) differ 
considerably from the idealized "ionic" configurations, 
and the respective deviations increase with the d-p 
covalent overlap in the oxides. 

IV. CONCLUSIONS 

In summary, we have performed FLAPW calcula- 
tions in order to systematically study ground state 
properties and electronic structure features in 3d and 
4d transition-metal oxides SrM03, where M = Ti, V, 
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Zr and Nb. The energy gap between 2p oxygen and d 
metal bands increases with the decrease in the covalent 
p-d overlap. The optimized lattice parameters, bulk 
modules, densities of states, band structures and charge 
density distributions are obtained and compared with 
the available theoretical and experimental data. The 
analysis of the d band population reveals that the 
deviation of the d metal electron configurations from 



the formal ionic state (d° for SrTiC>3, SrZrC>3 and d 1 for 
SrVC>3 and SrNbOs) increases with p-d covalent overlap. 
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FIG. 2: Total and projected DOSs for the cubic perovskites SrB0 3 (B = Ti, V, Zr and Nb). 




FIG. 3: Charge densities maps for SrTiOa and SrVC>3 in (100) plane. 



